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Hypokalaemic periodic paralysis is typically associated with mutations of voltage sensor residues in calcium or sodium channels 
of skeletal muscle. To date, causative sodium channel mutations have been studied only for the two outermost arginine residues 
in S4 voltage sensor segments of domains I to III. These mutations produce depolarization of skeletal muscle fibres in response 
to reduced extracellular potassium, owing to an inward cation-selective gating pore current activated by hyperpolarization. Here, 
we describe mutations of the third arginine, R3, in the domain III voltage sensor i.e. an R1135H mutation which was found in 
two patients in separate families and a novel R1135C mutation identified in a third patient in another family. Muscle fibres from 
a patient harbouring the R1135H mutation showed increased depolarization tendency at normal and reduced extracellular 
potassium compatible with the diagnosis. Additionally, amplitude and rise time of action potentials were reduced compared 
with controls, even for holding potentials at which all Nav1.4 are fully recovered from inactivation. These findings may be 
because of an outward omega current activated at positive potentials. Expression of R1135H/C in mammalian cells indicates 
further gating defects that include significantly enhanced entry into inactivation and prolonged recovery that may additionally 
contribute to action potential inhibition at the physiological resting potential. After S4 immobilization in the outward position, 
mutant channels produce an inward omega current that most likely depolarizes the resting potential and produces the hypo- 
kalaemia-induced weakness. Gating current recordings reveal that mutations at R3 inhibit S4 deactivation before recovery, and 
molecular dynamics simulations suggest that this defect is caused by disrupted interactions of domain III S2 countercharges with 
S4 arginines R2 to R4 during repolarization of the membrane. This work reveals a novel mechanism of disrupted S4 translocation 
for hypokalaemic periodic paralysis mutations at arginine residues located below the gating pore constriction of the voltage 
sensor module. 
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Hypokalaemic periodic paralysis mutations disrupt S4 translocation 

Introduction 

Hypokalaemic periodic paralysis is an autosomal dominant skeletal 
muscle disorder named for recurrent episodes of flaccid limb weak- 
ness with reduced ictal serum potassium. Triggering factors for these 
episodes include carbohydrate-rich meals, rest after exercise, early 
morning hours, and emotional stress (reviewed by Fontaine, 2008; 
Cavel-Greant a/., 2012). Nearly all causative mutations in hypo- 
kalaemic periodic paralysis are associated with the Cav1.1 calcium 
channel encoded by CACNA1S or the Nav1.4 sodium channel 
encoded by SCN4A (reviewed by Matthews et al., 2009). For 
either of these channels, mutations are located at one of two argin- 
ines (R1 or R2) nearest the extracellular end of voltage-sensitive S4 
segments, whose movement in response to membrane depolariza- 
tion is linked to opening of the central conducting pore. Functional 
expression of channels with R1 or R2 mutations reveals an aberrant 
leak current through the S4 gating pore (Sokolov etal., 2007, 2010; 
Struyk and Cannon, 2007; Struyk et al., 2008; Jurkat-Rott et al., 
2009; Francis et al., 2011). This so-called 'omega' current allows 
cation influx at the resting state of the sodium channel to depolarize 
the muscle membrane. Additional depolarization, with reduced 
conductance of the inwardly rectifying potassium channels by hypo- 
kalaemia, impairs the generation of action potentials and leads to 
muscle weakness (Cannon, 2010; Catterall, 2010; Tricarico and 
Camerino, 2011; George etal., 2012; Jurkat-Rott etal., 2012). 

Recently, a sporadic case of hypokalaemic periodic paralysis 
with Nav1.4 mutation R1135H in the domain III voltage sensor 
was described (Sung et al., 2012). Until now, functional expres- 
sion has not been published for this mutation. However, based on 
its location (R3) and hypothesized function of the affected S4 
segment, the mutant residue should traverse the gating pore of 
the membrane where the omega current is thought to originate, 
but only at depolarized potentials in activated or inactivated states 
of the channel (Sokolov et al., 2005). Several R3:DIIS4 mutations 
in normokalemic periodic paralysis cause an outward omega cur- 
rent at depolarized voltages, which is sustained with inactivation 
of the channel (Sokolov et al., 2008). In this work, we describe 
R3:DIIIS4 Nav1.4 sodium channel mutations found in hypokal- 
aemic periodic paralysis that promote inward and outward 
omega currents under different conditions. The clinical presenta- 
tions for two families with R1135H, and one with a novel R1135C 
mutation, are compared to the biophysical defects caused by these 
mutations as observed with electrophysiological recordings from 
patient muscle fibres or in heterologous expression, and as 
simulated with molecular dynamics trajectories. 



Materials and methods 

Families 

Whole EDTA blood was taken from patients for CACNA1S and SCN4A 
analyses. Complete gene sequencing was performed using PCR ampli- 
fication and Sanger sequencing as described previously (Jurkat-Rott 
et al., 1994, 2000). Patient studies were approved by the institutional 
review board in Dim and conducted according to the Declaration of 
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Helsinki. Informed consent was obtained from patients, and from 
volunteers with no evidence or history of muscle disease. 

Magnetic resonance imaging 

^H- and ^^Na-AARI of the lower legs of a patient from the HypoPP317 
(R1135H) family were performed as described previously (Jurkat-Rott 
et al., 2009). 

In vitro studies on native muscle fibres 

Muscle specimens were removed from the quadriceps muscles of the 
index patient of Family HypoPP4 (R1135H) and three adult individuals 
with no neuromuscular disease under regional anaesthesia. Muscle 
specimens were ~28mm in length and 5 mm in diameter. They 
were further prepared into small bundles and allowed to reseal over 
2 h in a solution also used for resting membrane and action potential 
measurements. This contained (in mM): NaCI 108, KCI 4, CaCb 1.5, 
MgS04 0.7, NaHCOB 26.2, NaH2P04 1.7, Na-gluconate 9.6, glucose 
5.5, sucrose 7.6, 290 mOsmol/l, maintained at 37°C, with pH adjusted 
to 7.4 by bubbling with 95% O2 and 5% CO2. In some experiments, 
extracellular was decreased to 1 mM. 

Membrane potentials were measured by use of a voltage amplifier 
(Turbo TEC-05, NPI Electronic Instruments) and glass microelectrodes 
filled with 3 M KCI and input resistances of 5-10 ^AQ. Histograms of 
the potentials were smoothed by density estimation (WARPing). The 
potentials exhibited a two or three-peak distribution of polarized and 
depolarized fibres and were displayed as probability density. 
Parameters were obtained from a Gaussian fit for each peak according 
to f(x) = yo + (A / w*sqrt(7r / 2)) x exp(-2 x ((x - Xn) / w) a 2) 
where yo is the minimum asymptote, A is density, and Xn is midpoint 
voltage. 

Action potential recordings were performed with a second micro- 
electrode. This electrode delivered a constant current to hold various 
resting potentials for at least 1 min to ensure recovery of voltage-gated 
sodium channels. Then, action potentials were elicited by short depo- 
larizing pulses. The maximum rate of rise was determined at the 
potential for which the slope of the upstroke of the action potential 
was maximal. These values were taken as index for the available 
sodium conductance (Cohen al., 1984) and used to determine 
slow inactivation parameters by Boltzmann fits. Data are presented 
as mean ± standard deviation (SD). Statistical significance was 
assessed by Student's f-tests for normally distributed data with a 
criterion of P < 0.05. 

Whole cell recording experiments 

Site-directed mutagenesis was performed using an overlapping PCR- 
based technique. Subsequently, the mutants were reassembled in the 
pRC/CMV plasmid (Invitrogen) for transfection by calcium phosphate 
precipitation in tsA201, a mammalian cell line. Standard whole cell 
recording methods were used, as previously described (Kuzmenkin 
et al., 2002). Capacity transients were eliminated by a -p/4 protocol. 
Series resistance errors were <3mV and corrections were made for 
liquid junction potentials. Data were filtered at 3 kHz, and acquired 
using pCLAMP 10 (Axon Instruments). Patch electrodes contained 
(in mM): CsF 105, NaCI 35, EGTA 10, Cs-HEPES 10. The bath 
contained NaCI 150, KCI 2, CaCb 1.5, MgCb 1, Cs-HEPES 10. 
Solutions were adjusted to pH 7.4 using CsOH, and experiments 
were performed at room temperature (20 to 22°C). 

Whole cell data for activation, inactivation, and recovery were ana- 
lysed by a combination of pCLAMP and ORIGIN (MicroCal). Data are 
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presented as means ± SEAA. Statistical significance was assessed by 
Student's ^-tests for normally distributed data with a criterion of 
P < 0.05. 

Cut-open oocyte experiments 

The rat isoform of SCN4A in pGEMHE was generously provided by Dr 
Steven Cannon (University of Southwest Texas Medical Centre, USA). 
Constructs rR1128H and rR1128C were made using Quik Change XL 
II site-directed mutagenesis kits (Agilent Technologies) and confirmed 
by sequence analysis of the coding region. Beta subunit SCN1B was 
housed in pgH19. Plasmids were linearized with Nhel {SCN4A) or 
Hindlll (SCNW), and messenger RNA was generated using T7 or T3 
RNA polymerase (mAAESSAGE mAAACHINE kits, Ambion). Transcripts 
were co-injected at a mass ratio of ^a : 3Pi into oocytes from adult 
Xenopus laevis frogs. Surgical procedures for extraction of oocytes 
were performed according to lACUC guidelines of the Animal Use 
and Care Committee at ISU. Oocytes were cultured in medium 
containing in mAA: NaCI 96, KCI 2, CaCb 1.8, AAgCb 1, HEPES 5 
and Na pyruvate 2.5, with lOOmg/l gentamicin sulphate and 4% 
horse serum (Hyclone Laboratories, Fisher Thermo Scientific) at pH 
7.4, for 2 to 6 days before recording. All salts and antibiotics were 
obtained from Sigma Chemical Co. 

To record gating and leak currents, we used a CA1B amplifier 
(Dagan Corporation) with PULSE 8.67 software (HEKA). Oocytes 
were placed between top and guard chambers. The external solution 
contained in mAA: NAADG (A/-methyl-D-glucamine) 115, Ca(OH)2 1.5, 
BaOH 2.5 and HEPES 10, pH 7.4 adjusted with methanesulphonic acid 
(Sigma). For external ionic substitution, NAADG was replaced by the 
hydroxide of Na"^, Li ^ or 115, or partially by guanidine sulphate 
60. The internal solution contained in mAA: NAADG 115, EGTA 10 and 
HEPES 10, pH 7.4, with ionic substitutions the same as above. Before 
recording, the oocyte was permeabilized with 0.1 to 0.5% saponin to 
control the electrical and ionic interior of the oocyte. Experiments were 
started 40 to 50min later, and ionic currents were blocked with 2|iAA 
tetrodotoxin. 

Raw leak currents were recorded with R3 of IIIS4 below or above 
the gating pore constriction site. To separate background leak from 
omega currents, the mean current amplitude between 60 to 95 ms of 
the test potential was taken, and the linear component of the leak 
current was subtracted from the original l/V relation. For each experi- 
ment, leak current amplitudes were normalized against anodic gating 
charge displaced in response to a 20 ms step depolarization to 
+ 40mV with leak subtraction using -p/8. 

In a separate set of experiments the anodic gating charge associated 
with activation (IgoN) leading to fast inactivation, and cathodic gating 
charge associated with subsequent deactivation (Igopp) were mea- 
sured. The integrals of the gating currents IgoN and Igopp were deter- 
mined and expressed as Igopp / IgoN, as a measure of the return of 
the available gating charge during S4 translocation preceding recovery. 
The extent of gating charge return during hyperpolarization was quan- 
tified as mean Igopp / IgoN as a function of command voltage. 

Molecular dynamics simulations 

A homology model of the rNav1.4 domain III voltage sensor module 
(S1-S4) was constructed using the A chain of the bacterial sodium 
channel NavAb as template (3RVY.pdb, Payandeh et al., 2011), in 
AAODELLER, as described in Groome and Winston (2013). The 
rNav1.4 voltage sensor module was embedded in a POPC (1-palmitoyl 
2-oleoyl-sn-glycero-3-phosphocholine) lipid bilayer bounded by aque- 
ous layers of 200 mAA NaCI using VAAD (visual molecular dynamics). 



and NAAAD (not another molecular dynamics) was used for molecular 
dynamics simulations. Initially, all components of the model except for 
lipids were constrained for 0.5 ns. The system was then minimized for 
another 0.5 ns, followed by an equilibration step with protein con- 
strained before 10 ns of free equilibration at a voltage bias of 
-1200mV, similar to conditions described in Delemotte et al. 

(2011) , Gosselin-Badaroudine et al. (2012) and Khalili-Araghi et al. 

(2012) . Selected frames were saved as .pdb files. Interatomic distances 
between DIIIS2 negative countercharges D1066 or El 076 and DIIIS4 
arginine residues R2, R3 or R4 were measured with Jmol 
(Supplementary Table 1). 

Results 

Families with hypokalaemic periodic 
paralysis mutations 

HypoPP4 family with heterozygous Nav1.4 
R1135H (C.C3404A) mutation 

The 47-year-old male proband had spontaneous paralytic attacks 
of up to 3-h duration, typically associated with slight hypokal- 
aemia of 3.3 mM K^. Onset was at age 16, and frequency was 
bi-monthly. Weakness was triggered by heavy meals, alcohol and 
stress. Oral K"^ salts ameliorated, and acetazolamide at 250mg/d 
prevented these episodes. Oral glucose provocation with 2g/kg 
together with 20 units of subcutaneous administration of insulin 
resulted in decreased serum K"^ from 4.1 to 3.4 mM (normal range 
3.5 to 5.0 mM) and generalized muscle weakness lasting for 
several hours. In early adulthood, a fluctuating weakness of the 
upper limbs developed. Frequency and severity of paralytic attacks 
decreased after age 40 when moving to the warmer climate of 
Polynesia. A muscle biopsy revealed vacuolar myopathy containing 
tubular aggregates and focal, type 2 fibre atrophy. An excised 
specimen of muscle fibres was used for membrane potential meas- 
urements. Molecular genetics identified a heterozygous Nav1.4 
R1135H mutation. The mutation was found in his similarly 
affected mother, whereas it was excluded in each of two non- 
affected children. 

HypoPP317 family with heterozygous Nav1.4 
R1135H (C.C3404A) mutation 

At age 16 the 21 -year-old proband with unaffected parents had 
cortisol-induced paraesthesias and muscle stiffness, and weakness 
in the extremities followed by weakness of trunk and respiratory 
muscles. In the intensive care unit, serum K"^ was markedly 
reduced to 1.9mM and creatine kinase was slightly increased to 
386 U/l (normal < 270). Infusion of 100 mEq potassium in 2.51 
solution administered over a 24-h period completely resolved the 
weakness. 

In addition to endogenous spikes in Cortisol levels (circadian 
peak in early morning, and after stress) and iatrogenic Cortisol 
spike, other triggers were salt and liquorice, whereas cold 
temperature and exercise were not. In the 2 years after onset, 
episodes were bi-weekly and lasted 3 to 4h. Two episodes were 
clinically similar to the Guillain-Barre syndrome and required hospi- 
talization and intravenous administration because of 
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hypokalaemia of 2.3 and 3.0 mM, respectively. Since the last hos- 
pitalization, the patient has been taking 50mg/d epierenone in 
addition to daily potassium supplementation of 240 mEq. 
Molecular genetics identified a heterozygous Nav1.4 R1135H 
mutation that seemed to be de novo as it was excluded in his 
non-affected parents. 

''H-MRI and ^^Na-MRI imaging have been used to detect 
muscular degeneration, oedema and sodium overload in hypo- 
kalaemic periodic paralysis patients harbouring R1 or R2 mutations 
of sodium or calcium channels (Jurkat-Rott et al., 2009). In the 
present study, ''H-MRI of the R1135H patient at age 21 showed 
neither fatty degeneration nor muscular oedema. The ^^Na-MRI 
revealed a muscular sodium signal intensity of 0.97 when normal- 
ized to a 0.3% saline reference solution. This ratio was in the 
range of controls (/i = 10; mean value ± SD: 0.99 ±0.12, 95% 
confidence interval of 0.93-1.05) suggesting a weak effect of 
the mutation (Weber et al., 2006). 

HypoPP181 family with homozygous Nav1.4 
R1135C (C.C3403T) mutation 

The female proband had consanguineous healthy parents (first 
cousins). Her siblings had no muscle symptoms, and her daughter 
was not affected. The patient had the first episode at age 14, 
comprising a paralysis of neck and limb muscles, and respiratory 
difficulty with sudden onset and lasting 1 h. Initially the weakness 
episodes had a frequency of about three times per week, without 
having a clear precipitating factor. At a later stage in life, admin- 
istration of catapresan for hypertension worsened her muscle 
symptoms, and convertan elicited attacks of respiratory weakness. 
Serum was checked several times with values mostly in the 
normal range, but at least for two episodes levels of 2.9 and 
3.4 mM were noted, the latter despite a preceding oral K"^ inges- 
tion. In the EMG there was no myotonia, and motor and sensory 
nerve conduction velocities were normal. 

With time, episodic weakness became more frequent and some 
permanent weakness became evident. In the last years of her life 
she was wheelchair-bound most of the time. The permanent 
weakness motivated a biopsy of the biceps brachii muscle at age 
60, revealing myopathy without biochemical alteration. Molecular 
genetics identified a homozygous Nav1.4 R1135C mutation that 
was heterozygous in her non-affected daughter. Parents were 
unavailable for testing. 

Intracellular action potentials 

Action potentials in fibres harbouring the R1135H mutation 
revealed undershoots of - 10.8 ± 9.5 mV instead of overshoots 
as in control fibres (+10.5 ± 7.9 mV, P < 0.05; Fig. 1A). 
The maximum rates of rise of R1135H action potentials were 
~50% less than those from control fibres (Fig. IB). The 
normalized maximum rates of rise, i.e. the steady-state slow inacti- 
vation curves, showed similar voltage dependence for R1135H 
and control fibres. Boltzmann fits to these curves yielded param- 
eters of Vh slowo.5 = -71.0mV and slope factor k = 8.0 mV (wild- 
type), compared with values of -74.8mV and 9.3 mV 
for R1135H. 
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Resting membrane potentials 

In 4mM extracellular K"^, control and R1135H patient 
fibres showed a bimodal distribution of the resting potential. 
Gaussian fitting revealed peak values of PI = -81.1 ±0.3 and 
P2 = -58.7 ± 0.4 mV, similar to control fibres of 
PI = -80.9 ±0.1 and P2 = -60.2±0.4mV (Fig. 1C). However, 
more fibres were depolarized to P2 in the patient fibres as shown 
by the respective densities for PI and P2, which were 0.38 ± 0.03 
and 0.63 ± 0.05 compared with control fibres with values of 
0.80 ±0.03 and 0.04 ± 0.01 (Fig. 1C). 

Decreasing the extracellular K"^ to 1 mM led to a tri-modal 
distribution for both patient and control fibres. Under these 
conditions, the peak values were tendentially shifted to more 
depolarized potentials for the patient fibres, i.e. peak values of 
the three populations for patient fibres were PI = -86.4 ±0.8, 
P2 = -66.7 ± 0.5 and P3 = -51.1 ± 0.8, compared with the con- 
trol fibres with values of PI = -89.7 ± 0.6, P2 = -72.7 ± 0.2 and 
P3 = -52.2 ± 0.5 (Fig. ID). The density of the three peaks were 
0.03 ±0.02, 0.61 ±0.05 and 0.29 ± 0.05 for patients, and 
0.14 ±0.03, 0.68 ±0.03 and 0.25 ± 0.04 for controls, 
respectively. This can be interpreted as a borderline pathological 
result compatible with the diagnosis of hypokalaemic periodic par- 
alysis. The ability of patients' fibres to repolarize after strong 
depolarizations in chloride-free solution was not markedly 
impaired with PI = - 85.7 ± 0.4 mV versus - 79.9 ± 0.4 mV, 
P2=-65.1±1mV versus -63.7 ± 1.1 mV for control and 
patients, respectively (Supplementary Fig. 1). This is in contrast 
to the single peak at -50mV that fibres of hypokalaemic periodic 
paralysis patients show (Rudel etal., 1984). Taken together, these 
results suggest a rather weak effect of the mutation. 

Whole cell currents 

Figure 2A displays typical recordings of depolarization-induced 
sodium currents elicited in tsA201 cells. The corresponding volt- 
age-conductance curves according to a Boltzmann function are 
shown in Fig. 2B. Hypokalaemic periodic paralysis mutants dis- 
played reduced voltage sensitivity with k = -11.38 ± 0.28 
(R1135H) and -9.12 ±0.21 (R1135C) versus -7.54 ±0.16 
(wild-type) as expected for neutralization of a positive charge in 
the S4 voltage sensor. Additionally, activation of R1135C was 
significantly left-shifted by almost 4mV. In contrast, R1135H did 
not affect V0.5, but slowed activation kinetics by a factor of ~1.5 
(Fig. 2C). Activation is also slowed with histidine mutations in S4 
of domain II (Kuzmenkin et al., 2002), perhaps as a result of the 
bulky side group added by that substitution. 

Figure 2B shows steady-state inactivation curves for wild-type 
and mutant channels. Fits to a Boltzmann function revealed a 
significant left shift for both mutants, V0.5 = -97.9 ± 1 .36 mV 
(R1135H) and - 100.8 ± 1 .46 mV (R1135C) versus 
-84.7±1.07mV (wild-type), as well as a reduced slope 
k = 7.32 ±0.19 (R1135H) and 8.83 ± 0.43 (R1135C) versus 
5.80 ±0.18 (wild-type). Additionally, the extent of closed state 
inactivation was markedly increased for both mutants (Fig. 2D). 
Both of these phenomena are consistent with the observed slow 
recovery from the fast inactivated state (Fig. 2E) fit by a 
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Figure 1 (A) Representative action potentials elicited from R1135H patient or wild-type control muscle fibres from 1 min holding 
potential of -120mV. (B) Maximum rate of rise for R1135H patient {n = 6) and wild-type control {n = 6) fibres from varied holding 
potentials. (C) Bimodal distribution of membrane potentials of R1 135H patient fibres (n = 106) and wild-type control fibres (n = 256) with 
4 mM in the recording solution. Note the increased fraction of depolarized fibres in the patient. (D) Tri-modal distribution of membrane 
potentials of R1135H patient fibres {n = 55) and wild-type control fibres (n = 61) with decreased K"^ to 1 mM. Note the borderline 
depolarizing shift of the peak potentials in the patient fibres compared with controls. 



two-exponential function, whereby the faster component differed 
significantly: tfast = 7.47 ± 0.72 ms (R1 1 35H) and 9.57 ± 1 .26 ms 
(R1135C) versus 2.48 ± 0.22 ms (wild-type). Fig. 2F gives an 
overview of the time constants associated with fast inactivation 
for wild-type and mutant channels. 

Results from experiments testing the effects of mutations on 
slow inactivation are shown in Supplementary Fig. 2. In wild- 
type channels, the midpoint of the steady-state slow inactivation 
curve Vsi 0.5 was -68.5 ± 1 .38 mV, with a slope factor of 
12.3±0.7mV (Supplementary Fig. 2A). R1135H produced a 
4mV left shift of Vsi 0.5 (-72.2 ± 1 .3 mV, slope factor of 
13.9 ± 0.98 mV), and R1135C produced a 6mV left shift of Vs, 
0.5 (-74.6 ± 1.46 mV, slope factor of 15.4 ± 1.11 mV). Wild-type 
channels entered into slow inactivation with a time constant of 
4.37 ± 0.29s, with 87% completion, compared with R1135H 
(4.44 ± 0.39 s, 91% completion) and R1135C (3.56 ± 0.31s, 
92% completion; Supplementary Fig. 2B). Recovery from slow 
inactivation was prolonged in hypokalaemic periodic paralysis 
mutants (Supplementary Fig. 2C), with time constants of 



3.60 ± 0.52 s (R1135H) or 3.18 ± 0.35 s (R1135C), compared 
with 1.36 ± 0.08 s in wild-type channels. 

Cut-open oocyte recordings 

First, we measured leak currents in wild-type and mutant channels 
in experiments that compared the effect of external or internal 
substitution of small cations, guanidinium, or the large cation 
NMDG on outwardly or inwardly directed currents over the volt- 
age range of -160mV to +40mV, from a holding potential of 
-lOOmV. With Na^, K"^ or Li"^ ion in the external solution, we 
did not observe an inwardly directed omega current similar to that 
for hypokalaemic periodic paralysis mutations at R1 or R2 in do- 
mains I to III of Navl .4 (data not shown). To test for an outwardly 
directed current in mutant channels, we performed experiments 
with cationic internal substitution (Fig. 3). rR1128H/C promoted a 
cation-selective, outwardly directed gating pore current with in- 
ternal substitution of 60 mM guanidinium, and smaller outward 
currents with 115mM Na"^, K"^ or Li"^ ion. These currents were 
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Figure 2 (A) Sodium currents elicited by depolarization to voltages ranging from -85 to +55 mV from a holding potential of -140mV 
as recorded from tsA-201 cells expressing wild-type hNav1.4, R1135H, or R1135C channels. Calibration: horizontally, 5 ms; vertically, 
wild-type hNav1.4, R1135C 2 nA, R1135H 1.5 nA. (B) Conductance-voltage relationships from Boltzmann fits of steady-state activation 
and inactivation {n = 9-12). Steady-state fast inactivation was tested by a series of 300 ms depolarizing prepulses from -160 to -30mV 
followed by a test pulse to -lOmV. (C) Kinetics of activation determined by the 10-90% rise time of sodium currents between -60 and 
+ 30 mV. (D) Time course of closed state inactivation for durations from 0.05 to 300 ms at - 90 mV from a holding potential of - 1 40 mV 
whereby a test pulse to -lOmV determined the fraction of non-inactivated channels. (E) Recovery from fast inactivation from a holding 
potential of - 120 mV after an inactivating prepulse at - 10 mV for 100 ms. (F) Voltage dependence of time constants of fast inactivation: 
for -140 to -lOOmV, recovery from fast inactivation; for -100 to -70 mV, closed-state inactivation; for > -70 mV, open-state 
inactivation. 



larger for rR1128C than for rR1128H and showed rectification at 
voltages more positive than -60mV. With internal substitution of 
NMDG, or with oocytes injected with wild-type rNav1.4, leak 
currents exhibited a linear voltage dependence. 

These experiments suggested that R3 is located below the 
gating pore in the resting state of the channel. We wished to 
test whether R3 is above the gating pore in the activated/inacti- 
vated states of the channel, to investigate the potential role of an 
omega current during recovery from fast inactivation for hypokal- 
aemic periodic paralysis mutations. To move R3:DIIIS4 from below 
to above the gating pore constriction, we preconditioned oocytes 
to a holding potential of +40mV for 100 ms before testing for 
leak currents over a voltage range of -160mV to +40mV 
(Fig. 4). With external substitution of 115mM Na^ or 60 mM 
guanidinium, channels with rR1128H/C mutations produced 



apparently inwardly rectifying currents, with the cysteine mutation 
producing larger currents. 

In Fig. 4B, it is difficult to observe outward omega current be- 
cause gradients are reversed compared with Fig. 3. To clarify 
whether the conditioning prepulse additionally occluded the 
omega pore, we examined raw leak current amplitudes with ex- 
ternal sodium or guanidinium versus NMDG at -20mV to 
+ 40mV for wild-type and mutant channels without normalization. 
In wild-type rNav1.4 channels, leak currents at depolarized volt- 
ages were increased by 1.35-fold (sodium) and 1.52-fold (guani- 
dinium), compared with currents elicited in NMDG. Relative leak 
current increases in sodium or guanidinium of 1.85- and 3.49-fold 
(rR1128H) and 2.57- and 3.05-fold (rR1128C), respectively, sug- 
gest that there is some population of mutant channels following 
preconditioning that have an open omega gating pore. At more 
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Figure 3 (A) Cut-open oocyte recordings of wild-type rNav1.4 and hypokalaemic periodic paralysis mutant channels with internal 
substitution of 115 nM NMDG. Traces are every other sweep showing non-subtracted currents in response to command potentials from 
-160mV to +40mV, from a holding potential of -lOOmV. (B) Responses to same protocol, after internal substitution with 60 mM 
guanidinium / 60 mM NMDG. (C) Plots of normalized leak currents for wild-type and mutant channels for internal substitution of NMDG, 
monovalent cations or guanidinium. Currents were normalized with respect to gating charge integral recorded in response to +40mV 
depolarization. 



hyperpolarized potentials, it is likely that additional R3:DIIIS4 
mutant residues enter that pore to produce the larger inwardly 
directed currents. We interpret these findings in the sense that 
the S4 movement during channel recovery is increasingly 
disturbed. To further examine this, we studied gating charge 
movement in wild-type and mutant channels, and modeled the 
molecular dynamics of the S4 translocation. 

On and off gating current recordings 

We performed gating current measurements of wild-type rNavl .4 
and R3C/H using a protocol similar to that used in experiments 
shown in Fig. 4, to directly measure S4 deactivation with hyper- 
polarization from the inactivated state. We did this to investigate 
that possibility that hypokalaemic periodic paralysis mutations at 
R3:DIIIS4 slow the deactivation following fast inactivation, to ex- 
plain the effect of these mutations to prolong recovery of 
channels. 

The integral of the gating currents (charge) (Fig. 5A) was com- 
pared as IgoFF / IgoN for hyperpolarizing commands of durations 
up to 10 ms. Figure 5B shows the time course of deactivation, with 
command to -105mV. An initial rapid rise in IgoFF / IgoN was 



observed for wild-type and mutant channels, reflecting the ~40% 
of gating charge that is unaffected (non-immobilized) by fast 
inactivation (Cha et al., 1999). The slower phase of deactivation 
was significantly prolonged in R3C/H, compared with wild-type 
channels. 

The extent of the inhibition of S4 translocation in mutant 
channels compared with wild-type rNavl .4 is shown in Fig. 5C, 
with measurements taken as mean IgoFF / IgoN fi'om 9.2 to 10 ms 
over the voltage range of -155 to -55mV. Significantly less 
charge was returned in the mutant channels compared to wild- 
type for hyperpolarizing commands more negative than -70mV; 
this effect produced by R3C/H is also observed in the current 
traces at 10 ms pulse duration shown in Fig. 5A. Importantly, 
our results show that hypokalaemic R3:DIIIS4 mutations disrupt 
the inward S4 translocation expected during recovery of channels 
from the fast inactivated state. 

Molecular dynamics simulations 

We used molecular dynamics simulations as an additional ap- 
proach to investigate translocation of the voltage sensor of 
domain III in wild-type rNavl .4 and hypokalaemic periodic 
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Voltage (mV) Voltage (mV) Voltage (mV) 

Figure 4 (A) Cut-open oocyte recordings of wild-type rNavl .4 and hypokalaemic periodic paralysis mutant channels during recovery from 
fast inactivation. From a holding potential of - 1 00 mV, the oocyte membrane was held at + 40 mV for 1 00 ms, before command test pulses 
at voltages from - 1 60 mV to + 40 mV. Traces shown (every other sweep) are responses in experiments with external solution of 1 1 5 mM 
NMDG, 115 mM Na^ , or 60 mM guanidinium. (B) Plots of normalized current amplitudes for each of the conditions shown in A. 



paralysis channels, with conditions similar to those for 
electrophysiological experiments shown in Fig. 4. We did this to 
investigate whether the electrostatic interactions of DIIIS4 positive 
charges with negative countercharges suggested by recent studies 
(Gosselin-Badaroudine et al., 2012; Groome and Winston, 2013) 
are disrupted by the hypokalaemic periodic paralysis mutants at 
R3:DIIIS4 in a manner consistent with the findings from our func- 
tional studies. Figure 6A shows locations of relevant residues in the 
homology model of S1-S4 for domain III of rNav1.4. After equili- 
brating the protein in a POPC membrane in silico (Fig. 6B), nega- 
tive potential was applied and the simulations were run for 10 ns 
(Fig. 6C). At the time points shown, proximities of side chains for 
DIIIS2 countercharges D1066 and E1076 with R2, R3/H3/C3, and 
R4 were measured, and are given in Supplementary Table 1. 

For wild-type rNavl .4, R2 and R3 residues in S4 were each in 
close proximity (<3A) to upper negative countercharge D1066 in 
S2 at the start of simulation, as was R4 with the lower negative 
countercharge El 076 (Fig. 6C, top). At 4 ns into simulation, R2 
and R3 were 6 A or more distant from the upper S2 counter- 
charge. Subsequently, the side chain of R3 passed the aromatic 
ring of F1073 to reside near El 076 in S2. At 10 ns, R4 was no 
longer in proximity to El 076. 

MD simulations for rR1128H/C shown in Fig. 6C (middle and 
bottom panels) revealed several differences in countercharge to 



DIIIS4 interactions, compared to wild-type channels. For example, 
R3C and R3H were each 6 A or more distant from D1066 at the 
start of simulation, and neither came in close proximity to El 076 
in response to membrane hyperpolarization. For each R3 mutant 
channel, R2:D1066 and R4:E1076 proximities were similar to 
those for wild-type rNavl .4, at the start of simulation. However, 
in rR1128C the R2:D1066 proximity measurement remained 
within 3 A throughout simulation, and for both mutants, R4 re- 
mained in close proximity to El 076. 

Discussion 

Genotype-phenotype correlation 

Our results with native R1135H muscle fibres show that mutation 
at R3 increases the probability of depolarized resting membrane 
potentials under normal conditions. At low external potassium and 
after strong depolarizations in chloride-free solutions, the depolar- 
ization tendency is only marginal. Likewise, ^^Na-MRI showed a 
normal result. Additionally, the amplitude of the inward omega 
currents were only ~5% to 10% to those of R1 or R2 previously 
described (Sokolov etal., 2010). Furthermore, since inward omega 
currents only occur when IIIS4 has moved outward and been 
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Figure 5 (A) Protocol used to measure on and off gating currents. Tetrodotoxin-blocked wild-type and mutant channels were depolarized 
for 30 ms to + 40 mV to elicit IgoNr followed by a series of hyperpolarizing commands with variable voltage ( - 55 mV to - 1 55 mV) and 
variable duration (0 to 10 ms in 0.2 ms intervals) to elicit IgoFF- Traces shown are for 10 ms command to -105mV. (B) Kinetics of 
deactivation at - 105 mV for pulse durations from 0 to 10 ms, expressed as the ratio of IgoFF / IgoN- (O Voltage dependence of the effect 
of R3C/H to inhibit deactivation, measured at 9.2 to 10 ms hyperpolarizing pulse durations. 



immobilized, they are not produced by the total mutant channel 
population but only by a portion of these, further reducing the 
effect of the mutation. Therefore, the effects of both mutations 
may indeed be slight. 

As to the recessiveness of R1135C, we hypothesize that epigen- 
etic factors or the differential effects of the two mutations on 
activation and recovery of the central pore current may be decisive 
for development of symptoms and the resulting mode of inherit- 
ance. As the dominant mutation R1135H slows all transitions pre- 
sumably because of the bulkiness of the histidine residue, it could 
have a larger effect on the action potential and excitability in 
general than the recessive mutation R1135C. 

Mechanisms for the effects of 
R3:DIIIS4 mutations 

Effect of the outward omega currents 

Recordings from patient muscle fibres showed an overall inexcitabil- 
ity, reduced amplitude and slowed rate of rise of action potentials. 
The reduced rate of rise to 50% action potential amplitude cannot 
be explained by depolarization-induced sodium channel inactivation, 
since we repolarized these fibres down to -120mV for at least 
1 min prior to action potential measurements. Cut-open oocyte 
recordings of R3:DIIIS4 mutations showed an outward omega 



current that is activated over the depolarized voltage range, a find- 
ing that is in agreement with earlier studies showing that mutations 
at R3:DIIS4 that cause normokalemic periodic paralysis (Vicart et al., 
2004) also produce an outwardly directed omega current (Sokolov 
et al., 2008). Our action potential recordings therefore suggest that 
outward omega current in R3C/H may contribute to reduced excit- 
ability by shunting excitatory current through the gating pore in the 
depolarized voltage range for which action potential initiation is 
expected. In addition to its effect on action potential initiation, 
the outward omega current in R3C/H may provide a counteraction 
to chloride loss, explaining the ability of these fibres to repolarize in 
chloride-free solutions. This current may augment the intrinsic repo- 
larizing mechanisms of ATPase and NKCC transporters (Lehmann- 
Horn and Jurkat-Rott, 1999). Possibly, the outward omega current 
is also able to reduce ictal hypokalaemia down to the mild levels 
found in these patients. 

Effect of inward omega currents and disrupted 
S4 translocation 

We found that each of the hypokalaemeic periodic paralysis mu- 
tations at R3:DIIIS4 impedes the downwards movement of S4, or 
deactivation, when channels were inactivated, and thus S4 is 
immobilized. The resulting small inward omega current is the 
basis for typical pathogenesis of hypokalaemic periodic paralysis 
as previously described (Cannon, 2010; Catterall, 2010; Tricarico 
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Figure 6 (A) Homology model of wild-type rNav1 .4 DIM (domain III) voltage sensor module highlighting S2 and S4 residues; the R3 locus 
is in light green. (B) VSM model incorporated into POPC membrane. (C) Select frames at times shown for simulation of response to 
membrane hyperpolarization. Residue R3 (white arrow) is highlighted as light green (wild-type rNav1.4, top), orange (rR1128H, middle) 
or magenta (rR1128C, bottom). At each time, distance measurements were made for R2:D1066, R3/H3/C3:D1066, R3/H3/C3:E1076, 
and R4:E1076. 



and Camerino, 2011; George, 2012; Jurkat-Rott et al., 2012). 
However, for the mutations at R3 the omega current is produced 
by a different mechanism than the previously described mutations, 
namely by a disturbed S4 movement during recovery from the 
inactivated state. Specifically, whereas S4 translocation and R3 
passage upwards through the gating pore is relatively favourable, 
the reverse translocation is impeded and results in a prolonged 
non-occlusion of the pore i.e. enhancing the inward omega cur- 
rent. The effect of the hypokalaemic periodic paralysis mutations 
on S4 translocation may be their predominant effect to decrease 
membrane excitability in channels that have fast inactivated in the 
process of firing action potentials. For example, it has been shown 
that the voltage sensors in domains III and IV are immobilized 
during fast inactivation (Cha et al., 1999), with immobilization 
of the DIIIS4 sensor coupled to binding of the IFMT particle to 
its receptor (Sheets et al., 2000). Therefore, our finding that hypo- 
kalaemic periodic paralysis mutations at R3:DIIIS4 inhibit S4 trans- 
location from a state of immobilization is consistent with the 
finding that these mutations slow the overall recovery of channels 
from fast inactivation. 



Molecular dynamics simulations support the results of gating 
current measurements that show R3C/H impedes S4 translocation 
during deactivation. Proximity measurements taken during these 
trajectories suggest that DIIIS2 countercharges may be an import- 
ant substrate of voltage sensor translocation during recovery, and 
that the effect of hypokalaemic periodic paralysis mutations in- 
volves a decreased interaction of those S2 countercharges with 
R3 as the S4 segment traverses the gating pore towards the rest- 
ing state. This hypothesis is supported by a previous study using 
molecular dynamics to investigate intermediate states of voltage 
sensor translocation in rNav1.4 (Gosselin-Badaroudine et al., 
2012). In that simulation of wild-type rNav1.4, interaction of 
R3:DIIIS4 with the outer S2 charge (D1066) is proposed for the 
activated state of the channel, whereas interaction of R3:DIIIS4 
with an inner S2 charge (El 076) is proposed for intermediate and 
resting states. Additionally, R4 interacts with El 076 in the acti- 
vated and intermediate, but not resting state. Each of these inter- 
actions was reiterated in our simulations of rNav1.4. By 
incorporating the hypokalaemic periodic paralysis mutations into 
our simulations, we found that R3H/C each may have limited 
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interaction with the outer and inner S2 residues during voltage 
sensor translocation, and that the R4:E1076 interaction persists 
in the resting state of the channel. Therefore, we propose that 
R3:DIIIS4 hypokalaemic periodic paralysis mutations disrupt native 
electrostatic interactions in the voltage sensor module necessary 
for recovery. 
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